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Abstract 

We investigate the pp — )• ppe'^e^ and quasifree pn pne^e^ reactions within an effective La- 
grangian model at the laboratory kinetic energy of 1.25 GeV for which experimental data have 
recently been reported by the HADES Collaboration. The model uses a meson-exchange approx- 
imation to describe the initial nucleon-nucleon (NN) scattering. Contributions to the reaction 
amplitudes are included from the NN bremsstrahlung as well as from the excitation, propagation 
and radiative decay of the A(1230) isobar state. It is found that the HADES data on the e'^e~ 
invariant mass distribution in the pp — > ppe~^e~ reaction are reproduced excellently by our model 
where the A isobar term dominates the spectrum. In the case of the quasifree pn — > pne~^e~ 
reaction, a strong sensitivity to the pion electromagnetic form factor is observed which helps to 
bring the calculated cross sections closer to the data in the higher dilepton mass region. 

PACS numbers: 25.75.Dw, 13.30.Ce, 12.40.Yx 
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Dileptons (e"'"e~) provide a valuable tool to investigate the properties of the strongly 
interacting matter at high temperature and density formed in the relativistic and ultra- 
relativistic nuclear collisions because after their production they travel to the detectors 
almost undisturbed by the surrounding baryonic matter. A recurring feature of the dilepton 
spectra measured in such collisions at low (DLS and HADES fl, 2]), intermediate [ super 
proton synchrotron (SPS) [3]] and high (PHENIX [^]) energies has been the significant 
enhancement observed in the intermediate dilepton mass region over the contributions from 
the electromagnetic decays of hadrons and long-lived mesons. While the major part of the 
excess yield seen at the SPS energies is attributed to the leptonic decay of the p meson 
(formed in the tt+tt^ annihilation process) with strongly modified spectral function in the 
dense and hot hadronic matter at relativistic heavy ion collider (RHIC) energies it is 
believed to be more due the strong thermal contribution from the partonic phase (see, e.g., 
Ref. B). 

At lower beam energies (1-2 GeV/nucleon), various transport models '7H10| have been 



unable to explain the DLS data which must be attributed to some inherent problems in the 
theory as the new measurements of the HADES Collaboration at these beam energies have 



confirmed the old DLS data 



Unlike the situation at high beam energies, the causes of 



this discrepancy - for the light systems at least - are unlikely to be related to the in-medium 
effects. The insufficiently known cross sections for the dilepton production in elementary 
proton-proton {pp) and proton-neutron {pn) collisions, are an important reason behind this. 
Indeed, in a recent transport model calculation 12|] it has been shown that if the input pn 
bremsstrahlung cross sections (which are calculated within the soft-photon approximation 



model 



13|). are scaled up in an adhoc manner by factors of 3-4, the observed dilepton 



yields of both DLS and HADES experiments at beam energies of 1-2 GeV/nucleon can be 
reproduced. 

However, the microscopic models of dilepton production in elementary A^A^ reactions 
differ in their predictions of the pn bremsstrahlung cross section. While the calculations 
performed within the effective Lagrangian models of Refs. 14j-ll7l| do not support the larger 

favor the enhanced cross sections 



pn bremsstrahlung yields, those of the model of Refs. |18| 
implemented in Ref. 12|. Therefore, to provide a reliable constraint for the dilepton yields in 
elementary reactions, the HADES Collaboration has very recently performed measurements 
for the dilepton production in not only the pp reaction but also in the quasifree pn reaction 
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at 1.25 GeV beam kinetic energy 19|]. The latter was measured by colliding a proton target 
with a deuteron beam of kinetic energy 1.25 GeV/nucleon and by detecting fast spectator 
protons from the deuteron breakup in a dedicated forward direction. In Ref. 19j the data 
on both pp and quasifree pn reactions were compared with the predictions of the model of 
Ref. [isl where it was noted that the calculations fail to describe the data for both reactions. 
While in the pp case the dilepton yields were overestimated in the entire range of the dilepton 
invariant mass (M), those of the quasifree pn reaction were overestimated (underestimated) 
at lower (higher) regions of M. 

The aim of this Rapid Communication is to investigate the dilepton production in pp 
and quasifree pn reactions at the beam energy of 1.25 GeV within the effective Lagrangian 



model (ELM) of Refs. Q, 17 1. In order to compare our calculations with the HADES 
dilepton yields, we have also considered the following additional features: (i) for the quasifree 
pn reaction the available energy in the center-of-mass (cm.) system has been smeared to 
include the momentum distribution of the neutron in the deuteron using the Argonne V18 
[2^ deuteron wave function. As a consequence, the dp reaction results in a smeared pn 
(quasifree) reaction where the available cm. energies could be in excess of the threshold 
for the rj meson production (see, e.g., Ref. 2l|), (ii) because of this we have included the rj 
Dalitz decay cross sections in the total theoretical yields for the pn reaction, and (iii) the 
contributions from the production and dileptonic decay of the subthreshold meson via the 
baryonic resonance A^*(1520), have been included for both pp and quasifree pn reactions. 

The Feynman diagrams [corresponding to both post emission and pre-emission (direct 
and exchange) processes] contributing to the dilepton production in the ELM are shown 
in Fig. 1. In case of the charged pion exchange (which happens for the pn reaction) the 
internal meson line can also lead to dilepton emission (Fig. Ic). Initial interaction between 
two incoming nucleons is modeled by an effective Lagrangian which is based on the exchange 
of the vr, p, u and a mesons. The coupling constants at the nucleon-nucleon-meson vertices 
are determined by directly fitting the T matrices of the A^A^ scattering in the relevant energy 



region 



These parameters are quite robust and have been used in successfu 



of NN NNn 



descriptions 



22], -> pAK+, pp pTPK+ |23|, l2J] and A^A^ NNrj ^ reactions. 
The dilepton production proceeds via excitation, propagation and radiative decay of 
the intermediate nucleon or resonance states at either of the two colliding nucleon vertices 
(Figs, la and lb). The nucleon intermediate states give rise to the A^A^ bremsstrahlung 
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FIG. 1: Graphical representation of our model to describe e~^e~ production in pp and pn collisions. 
Diagrams (a) and (b) represent the e^e~ emission from an external nucleon line, (c) e^e^ emission 
from an internal charged meson line in case of the pn reaction, (d) diagrams representing processes 
where e~^e~ is emitted directly from the charged meson-nucleon-nucleon-photon vertices. 



contribution. In calculations of various amplitudes we have used the same effective La- 
grangians (and the correspondin g p arameters) for all of the hadronic and electromagnetic 
vertices as those given in Refs. [la We recall that for the NNir vertex we have em- 
ployed a pseudoscalar (PS) coupling where no derivative term of the pion field is involved. 
Therefore, no extra term (corresponding to a contact or the seagull diagram) appears in the 
model at the NNn vertex when the electromagnetic coupling is included via the substitution 
—iemA'^, (m is +1, 0, -1 for positive, neutral and negative pions). However, at the 
NNp vertices such terms are always present as the corresponding coupling involve also the 



derivative of the meson field 



26 



28 



|. In any case, p-meson exchange terms contribute less 
than 5% to the total bremsstrahlung cross sections 29|. Thus, our calculations performed 
with a PS NNtc coupling, are almost free from the gauge invariance related ambiguities 
which could be associated with the contact terms (see, e.g., Ref. jsoj]). 

A note of caution should, however, be added here. With a PS NNn coupling the role 
of negative energy states may be overestimated, whereas these are quite suppressed with 
the corresponding pseudovector (PV) interaction. Nevertheless, only the bremsstrahlung 
dilepton production amplitudes are expected to be affected by the PS-PV coupling choice. 
Because, both pp — )■ ppe~^e~ and pn — )■ pne~^e~ reactions are dominated by the delta isobar 
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contributions at the beam energy of 1.25 GeV (as is shown in the following), our overall 
conclusions are not affected by the choice of the coupling at the NNn vertex. In any case, 
the nucleon-antinucleon vertex may be considerably suppressed compared to the NN vertex 



in presence of the timelike form factors 



3l|. 



As the beam kinetic energy of 1.25 GeV is below the threshold of the NNt] channel (1.258 
GeV), there is no contribution to the pp ppe^e~ cross sections from the r] Dalitz decay 
process. However, as explained earlier this can contribute to the quasifree pn reactions. 
The tt" Dalitz decay contributions must be taken into account as they dominate the cross 
sections at the lower ends of the dilepton invariant mass distributions in both the reactions. 

The 7] Dalitz decay is treated as a two step process - the rj meson production by reactions 
p + n p + n + rj and p + n — > d + rj, followed by the 77-meson Dalitz decay. The total 
cross sections for 77 meson production reactions have been taken from Ref. 25| where a good 
description of the corresponding experimental data [32| is obtained. The 77 Dalitz decay to 
'je~^e~ is calculated by using expressions given in Ref. [33|. A similar procedure is applied 
for the 77° case where the production cross sections have been taken from Ref. 22| while for 
its Dalitz decay the formulas of Ref. 33|] have been utilized. 

In calculations of dilepton yields from the production and decay of the subthreshold 
p° meson via the baryonic resonances, we consider only the A^*(1520) resonance - other 
higher lying resonances are expected to contribute negligibly at the beam energy con- 
sidered in this Rapid communication 3^. We suppose this reaction to proceed as a 
A^A^ — )■ -RA^ — i- p'^NN — e^e~NN process {R represents a resonance), which leads to 
the following factorization of the cross section. 



da(s,M) 
dM 



NN^NNe+e- 



da{s,M) 
dM 



NN-^p^NN 



.(M) 



P°*(M) 



where s is the square of the invariant mass associated with the incident channel. In Eq. (1) 
the first term on the right hand side represents the differential cross section for the meson 
production in A^A^ collisions which is calculated by following the procedure described in 
Refs. 35|, |36| using the same parameters as those given in Ref 36|]. The second term is the 
branching ratio for the p^ e^e~ decay. FpO^g+g- is the decay width of the p° meson to 
the dilepton channel, which is calculated within a strict vector meson dominance model as 
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in Ref. 3J]. r*°* is the total p-meson width which is given by [3 



r-(M) = ^P^^^^MjYTWy 

where A;^ = — m^. The parameter rc represents an interaction radius which is taken 

to be 2 fm. FpO^^Tr = 0.150 GeV. Eq. (2) represents the partial width for the p meson decay 
to the 27r channel only which has a branching ratio of nearly 100%. In our calculations, 
however, we have added to it also the width TpO^e+e-- 

In Fig. 2(a), we compare the calculated and the measured dilepton invariant mass (M) 
distributions for the pp — > ppe^e~ reaction at the beam kinetic energy of 1.25 GeV. We recall 
that for this reaction only diagrams 1(a) and 1(b) contribute. The theoretical cross sections 
have been folded with the appropriate detector acceptances provided to us by the HADES 



Collaboration [40| . We see that the total cross section (solid line) (obtained by the coherent 
summation of A^A^ bremsstrahlung and A isobar amplitudes which will be referred as QM) 
for this reaction is dominated by the A isobar terms (dashed line). The A^A^ bremsstrahlung 
contributions (dashed-dotted line) are smaller by almost an order of magnitude for lower 
values of M and by factors of 3-5 at higher M. The region of M < 0.15 GeV/c^ is dominated 
by the vr'' Dalitz decay cross sections. 

It clear that our QM cross sections for the pp — > ppe~^e~ reaction are in excellent agree- 



ment with the HADES data for M > 0.15 GeV/c^. In contrast to this, the model of Ref. 181 

n 

overestimates the data everywhere in this region (see Fig. 1 of Ref. [19!]). The similar ob- 
servation was also made in comparisons with the DLS pp dilepton data at 1.04 GeV beam 
energy in Ref. |l7|. Since this reaction is dominated by the A contributions, the larger A 
cross sections of Ref. [iS^ as compared to those of our model are the most likely reason for 
the differences seen in the predications of the two models. 

In Fig. 2(b), we show the same for the quasifree pn — )■ pne~^e~ reaction at 1.25 GeV beam 
kinetic energy. Because of the Fermi smearing, the tails of various contributions extend to 
M values larger than those of the pp reaction which is in agreement with the data. However, 
the shape of the pn spectra differs significantly in several ways from that of the pp case 
in the region of M beyond that dominated by the vr*^ Dalitz decay process. Firstly, the 
A^A^ bremsstrahlung contribution now is relatively larger although the A isobar term still 
dominates the total cross section. Secondly, the QM cross sections significantly underpredict 
the HADES quasifree pn data for M > 0.20 GeV/ which is in sharp contrast to the pp case. 
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FIG. 2: [color online] (a) The invariant mass distribution of the dileptons produced in the pp — >■ 
ppe~^e~ reaction at the beam energy of 1.25 GeV. The maximum allowed value of M is 0.545 GeV 
for this beam energy, (b) The same for the quasifree pn pne^e~ reaction at the same incident 
energy. The Fermi smearing as discussed in the text, has been applied to all the theoretical 
calculations shown in Fig. 2(b). Experimental data are from the Ref. 191 ]. 



The difference between theory and the data varies from factors of 2-3 at the lower mass values 
to more than an order of magnitude for M around 0.5 GeV/c^. It is important to understand 
this discrepancy between calculations and the data for the quasifree pn pne^e~ reaction 
as the dilepton excess in the intermediate mass range of 0.15 < MlGeV/c^) < 0.60 observed 
in the C + C collisions at 1 and 2 GeV/nucleon, can be explained by a superposition of 
experimental elementary pp and pn reactions as in-medium effects are almost negligible |l9|. 

In none of the results shown so far (in Figs. 2a and 2b as well as in Refs. 16|, Il7j ) 
electromagnetic form factors were considered at any of the vertices. However, in the earlier 
work reported in Ref. [l^ it was shown that the hadronic electromagnetic form factors 
have a significant influence on the dilepton spectra. Therefore, we now include the pion 
electromagnetic form factor (PEFF) [^^^(M^)] at the charged internal meson line (Fig. Ic). 
In the present exploratory study we use the same form factors at the pion and nucleon 
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vertices in order to preserve gauge invariance [2^ . We have used two parameterizations 
for F^(M^). The first one (to be referred as FFl) is written as 

2 

F (M'^) = ^ 

where rup is the p meson mass and Tp is the width for p — > vrvr decay. The assumption 
inherent in FFl is that the photon couples to the pion only via the p° meson. It reproduces 
the main features of the pion EFF both in time- and space- like regions (see, Ref. 38|]). The 
other, to be referred as FF2, is described extensively in Ref. 



^/,.9N 0.4 0.6 ml 

^-(^ ) = T7?7T? + 



MVA2 l-MV2m2m2-M2-imprp(M2)' 
where = 1.9 GeV^. The width rp(M2) appea ring in both FFl and FF2 has been calcu- 



lated by following the expressions given in Ref. [39|. FF2 is derived from the assumption 
that photon couples about 50% directly to the intrinsic quark structure of the pion and 
remaining 50% indirectly through the p° meson. FF2 provides a better description of the 
PEFF in the timelike region. The imaginary parts of both FFl and FF2 are proportional 
to two pion phase space - blow two pion production threshold both FFl and FF2 are real. 
It should be stressed here that we have put both form factors on the mass shell. Given 
the high virtuality of the internal pions, the form factors should be functions of both pion 
momentum and the momentum transfer. However, the knowledge about the off-shell pion 
form factor is still scanty - the recent extraction the PEFF from the JLab electroproduction 
experiments provide information about essentially the on shell pion form factors only (see, 
e.g., Ref. 4l|)- Therefore, we use the on shell PEFF in these calculations with a caveat that 
the off-shell PEFF could be larger than the on-shell one 42 1. 

In Fig 3(a), we show the effect of the PEFF and the contributions of the 7r° Dalitz decay 
and subthreshold p° decay process for the pp — > ppe^e~ reaction. We note that for this case, 
the introduction of the electromagnetic form factors (which are assumed to be the same for 
proton and pion vertices) makes hardly any difference to results obtained without them. 
The effect of FF2 type of PEFF is barely observed only at the extreme end of the spectrum. 
Results obtained with the FFl form factor are not shown here - there are even more closer 
to no PEFF results. Furthermore, the subthreshold p decay cross sections too are of some 
relevance only in the extreme tail region. 
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FIG. 3: [color online] (a) and (b) show the same reactions as in Figs. 2(a) and 2(b) but with elec- 
tromagnetic form factors included at pion and nucleon vertices. Also shown are the contributions 
of meson Dalitz decays and subthreshold decay processes. Total QM cross sections obtained 
with (FF2) and without (NEFF) electromagnetic form factors are shown by dashed and dashed- 
dotted lines, respectively. The simple sum of the meson Dalitz decays, decay and full quantum 
mechanical (with FF2 form factors) cross sections are shown by the full line. The data are taken 
from the Ref. Q- 

In Fig. 3(b), we show the total QM cross sections obtained without (NEFF) and with 
(FF2) electromagnetic form factor of FF2 type for the quasifree pn pne~^e~ reaction 
(where Fig. Ic also contributes) at the beam energy of 1.25 GeV. We have not shown 
explicitly the cross sections obtained with FFl type of PEFF in order not to overcrowd the 
figure - they lie between the NEFF and FF2 results. The larger cross sections obtained 
with FF2 form factors as compared to those with FFl can be traced back to the fact that 



38 



M- We note 



in the timelike region the former is significantly larger than the latter 
that with FF2 type of PEFF, the QM cross sections are significantly enhanced for M > 0.3 
GeV/c^ and are larger than t] and p° decay contribution by almost an order of magnitude at 
larger values of M. The t] Dalitz decay contributions drop off strongly for M beyond 0.50 
GeV/c^ due to phase space restrictions. In this region the p° decay cross sections become 
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relatively stronger. It is seen that the simple sum of the QM (with FF2) and the meson 
decay cross sections is able to reproduce the data now for M up to ^ 0.4 GeV/c^ and for 
M > 0.55 GeV/c^. It should, however be stressed that there is a danger of double counting 
by explicitly including p meson production and decay terms together with the form factor 
FF2 which implicitly includes a p meson bump. However, because the contributions of the 
explicit p° meson production process are relatively quite small as compared to that of the 
form factor FF2, this problem may not be too serious. 

We remark that the final state interaction (FSI) effects {pn and pp) estimated within 
the Watson-Migdal method increase the magnitudes of the cross sections with increasing M 
value. However, even at the extreme kinematical limits the FSI related enhancements in the 
cross sections are not more that 15-20% for these reactions at 1.25 GeV. This result is in 
agreement with those of Ref . 18| . Furthermore, the deuteronlike final states have also not 



been considered because the HADES measurements have ruled out such states in their data. 

In summary, we extended our effective Lagrangian model for dilepton production in 
nucleon-nucleon collisions by including the pion electromagnetic form factors at the internal 
meson line in a way that still preserves gauge invariance and employed it to describe the 
new data of the HADES Collaboration for these reactions. 

For the quasifree pn pne^e~ reaction, the inclusion of the electromagnetic form factors 
enhances significantly the cross sections for dilepton masses larger than 0.3 GeV/c^. This 
is for the first time that dilepton production data in elementary proton-neuron reactions is 
shown to be so sensitive to the pion electromagnetic form factors. Although this effect was 



already noted in the early work of Ref. ^J] but it could not be affirmed at that time because 
of the absence of data on the elementary pn process. We find that the simple sum of the 
tt" Dalitz decay and the ELM cross sections is able to describe the experimental invariant 
mass distribution of the dileptons everywhere except for the three points lying between 0.40 
-0.55 GeV/c^. The r] Dalitz and p° decay processes are of only minor consequence. 

For the pp — > ppe^e~ reaction, the ELM which remain almost unaffected by the elec- 
tromagnetic form factors, provides on its own a good description of the data for dilepton 
invariant mass > 0.15 GeV/c^. This is in a marked contrast to the results shown in Ref. lol 
where the model of Ref. |18| is found to grossly overestimate the data in this region. The sim- 
ple sum of ELM cross sections and those of the meson decay processes provides an excellent 
description of the data in the entire region of the dilepton invariant mass. 
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